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We propose a method to systematically control carrier densities at the interface of transition- 
metal oxide heterostructures without introducing disorders. By inserting non-polar layers sand- 
wiched by polar layers, continuous carrier doping into the interface can be realized. This method 
enables us to control the total carrier densities per unit cell systematically up to high values of 
the order unity. 
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1. Introduction 

Two-dimensional electron systems in semiconductor 
heterostructures have long been studied from scientific 
viewpoints as well as from industrial requirements, where 
carrier concentration can be controlled by the gate volt- 
age.^ In fact, physics of semiconductor developed in the 
past century has been scientific bases of modern electron- 
ics. In particular, two-dimensional interfacial microstruc- 
tures, such as metal oxide semiconductor field effect 
transistors (MOS-FET) and GaAs heterostructures, have 
played a major role. In these systems, fine-control tech- 
niques of two-dimensional electrons trapped at the in- 
terface have been developed. By using these structures, 
the sheet carrier density has been realized up to about 
^ lO^^cm"^ with remarkable suppression of disorder ef- 
fects. Quantum Hall effect^ has been observed in this re- 
gion of low carrier densities under strong magnetic fields. 

Meanwhile, clarification of correlation effects has been 
one of the most important issues in the condensed mat- 
ter physics. In the electron gas, the correlation effect 
becomes prominent for smaller density of electrons in 
general, because the averaged kinetic energy scales as 
rj^ for the mean distance between two neighboring 
electrons, while the mean Coulomb interaction energy 
scales as rj^. However, this rule does not necessarily 
hold when the electron density approaches a value of 
filling commensurate with the periodic potential of lat- 
tice formed from the atomic nuclei in the crystal. This 
commensurate filling can be reached, for example, at the 
density of one (or simple fractional number) per unit cell 
for conduction electrons, which is normally the order of 
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Fig. 1. (Color online) Schematic illustrations of reconstruction 
at interface. Dashed line plots (in red) represent the electronic 
potential along the [001] direction. Integers on the right of the 
layers represent charges of the layers, (a) The potential energy 
diverges along the [001] direction, (b) The polar instability can 
be avoided if Ti02 layer at the interface gets half an electron per 
unit cell. 



dimensional anisotropy such as transition-metal oxides, 
the on-site Coulomb repulsions strongly dominate phys- 
ical properties, when the sheet carrier densities reach 
10^'*-10^^cm~^ and the number of conduction electrons 
becomes close to a simple integer or fractional number 
per unit cell on average. Examples are found in the 
cuprate high-Tc superconductors, Mott insulators and 
anomalous metals.'^ The high carrier densities are real- 
ized in transition metal oxides conventionally by using 
chemical doping, which in most cases inevitably intro- 
duces disorder, unfortunately. Effects from disorder such 
as Anderson localization^ often obscure intrinsic corre- 
lation effects. On the other hand, one can realize carrier 
doping without disorder by FET.^ However, carrier den- 

13^ 



10^^ cm~3 in the conventional bulk crystal. This den- sity is achieved up to only about 10^-^cm ^ by using 



sity is much higher than the density of normal doped 
semiconductors and the above scaling naively suggests 
that the correlation effects are negligible in such ma- 
terials. However, the commensurability with the lattice 
leads to completely different physics, where strongly cor- 
related electrons emerge even at such high densities lead- 
ing to IVIott physics and charge ordering phenomena.'^ 
Indeed, in the strongly-correlated materials with two- 
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this method. If high carrier densities at transition-metal 
oxides could be realized without disorders, such systems 
could be ideal experimental stages of low-dimensional 
strongly-correlated electron systems and would make a 
major contribution to understanding of the correlation 
effects. 

Transition-metal oxide heterostructures are promis- 
ing candidates of realizing a disorder-free interface with 
higher density of carriers within an atomic scale thick- 
ness, for example by pulsed laser deposition.^ Recently, 
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Ohtomo and Hwang have reported the existence of a 
conducting electron layer at Ti02-LaO interface between 
two band insulators, SrTiOa and LaAlOa.® Other groups 
also reported the existence of superconductivity at the 
SrTiOg/LaAlOg interface.^ 

A mechanism for high conductivity at the transition- 
metal oxide interfaces has been proposed. The key 
idea is reconstruction of electronic distributions caused 
by a polar discontinuity, which is essentially understood 
from classical electromagnetism. Perovskite structures 
{ABO3) can be divided into alternating layers of AO 
and BO2 planes along the [001] direction. Sr^+0^~ and 
Ti^+Oj^ layers are charge-neutral layers, while in the 
ionic limit of LaAlOs, La'^+O^" has a positive charge and 
AP+02^ has a negative charge. A polar bilayer, such as 
La-'+O^^-AP+Oa^, makes no electric-field outside the 
bilayer, but gives a finite electric potential difference be- 
tween the plane above and below the bilayer. When the 
thickness of polar layers increases, the potential energy 
along the [001] direction increases, which yields an elec- 
trostatic instability (see Fig. 1 (a)). This instability can 
be avoided if 1/2 electron per unit cell is doped from 
the surface into the Ti sites at the interface (Fig. 1 (b)). 
These doped carriers contribute to electronic conduction. 
In this way, conductive interfaces made of transition- 
metal oxides can be realized, in principle, without dis- 
orders. 

This concept, the polar discontinuity, has been well 
appreciated since early times of studies on semiconduc- 
tors.^^' In Ge/GaAs interface, for example, Ge layers 
are charge-neutral, while GaAs layers have polarity along 
the [001] direction. To avoid the instability of the elec- 
tric potential from the polar layer of GaAs, half of the 
Ge sites at the interface are replaced with Ga (or As) 
atoms. Namely, the reconstruction occurs in the lattice 
systems on the atomic level. On the other hand, in the 
transition- metal oxide heterostructures, the reconstruc- 
tion is expected in the electron system. 

Here we note that, this mechanism enables us to dope 
only a discrete value of the carrier density at the inter- 
faces.^^' To investigate the nature of correlation effects, 
however, continuous tuning of the carrier densities is de- 
sired.'* 

In this paper, we propose a method to control the car- 
rier densities continuously without disorders by a specific 
procedure of doping into the transition-metal oxide het- 
erostructures. We demonstrate this charge controllability 
in the following chapter. 

2. Classical Calculation 

2.1 Doping into Interface 

In this chapter, we propose a method to systematically 
control carrier densities doped into the transition-metal 
oxide heterostructures without disorders. This method 
enables us to change the carrier densities continuously. 
Moreover, this method enables us to realize high densities 
of carrier doping. The idea of the systematic doping is 
very simple. We insert non-polar layers into the polar 
layers. Then, the doped carrier density changes, so that 
it is optimized to avoid the instability of the potential 
energy. In this paper, we refer to the non-polar perovskite 
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Fig. 2. (Color online) Schematic illustrations of reconstruction 
at interface with non-polar spacer. The (red and purple) dashed 
lines represent the electronic potential along the [001] direction. 
Integers on the right of the layers represent charges of the layers, 
(a) When a non-polar spacer is inserted into the polar layer, the 
potential difference between the interface and the surface does 
not change compared to the case of Fig. 1 (a) within the interval 
of the spacer, if there is no reconstruction, (b) On the other hand, 
in the case where half an electron is doped to the interface, the 
divergence of the electric potential is only incompletely avoided 
compared to the case of Fig. 1 (b). (c) Half doping into the 
interface is no longer the best way to avoid the divergence, and 
the optimum amount of the carrier filling changes from 1/2 to 
another value c. 

layers inserted into the polar layers as a spacer. 

An essence of the method can be understood in classi- 
cal physics. Let us consider non-polar perovskite, which 
has wide band gap and a lattice spacing close to the 
transition-metal oxide interface, as the spacer. Figure 2 
shows an interface where a spacer is inserted above a 
certain AIO2 layer. The electric potential along the [001] 
direction does not increase in the spacer region, if a re- 
construction does not occur (see Fig. 2 (a)). The poten- 
tial difference between the interface and the surface of 
the heterostructure is obviously the same in comparison 
with the non-inserted one (Fig. 1 (a)). However, in case 
that half an electron per unit cell is doped into the in- 
terface, the instability of the potential divergence along 
the [001] direction is only incompletely removed, because 
the electric potential increases in the region of the spacer 
(Fig. 2 (b)). In this case, doping 1/2 electron per unit cell 
into the interface is no longer the best way to avoid the 
potential divergence, and the optimum amount of carrier 
density changes from 1/2 to another value (Fig. 2 (c)). 

Now, we define L as the number of the polar layers, 
mB02 ^ the number of the spacer inserted above -BO2 
(e.g. Fig. 2 (a)), and m^o ^ the number of spacer in- 
serted above AO of the polar layer. We define a single 
layer (or a single spacer) as one unit cell along the [001] 
direction. The total energy per unit cell is obtained in 
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the classical electromagnetism as 

E = (iy+m/io+'7iB02)c'^ — (-^+2m^o)cx 1+const.. (1) 

where c is the doped carrier density per unit cell at the 
interface. In this section, we employ the energy unit by 
l/2e, where e is the dielectricity, for the sake of simplicity. 
Differentiating eq. (1) with respect to the carrier density 
c, we obtain the optimum number of the carrier density 
as 

1 L + 2mAo 

2 L + mAo + msos ' 
This result indicates that the carrier density can be con- 
trolled from to 1 by tuning the ratio of mB02 'tlao 
to L. Instead of the non-polar layers, if we can use ma- 
terials with other polarities as spacers, we obtain 



(2) 



{L + 2m,Ao)a + {~mAo + wbos)^ 



rriAO 



(3) 



mm I 2 ^ ^ 



where a is the polarized charge of the polar bilayer {+a, 
—a) and b is the charge of the spacer (+6, —6). Equation 
(2) represents the limit of the perfect polarity of the polar 
layer, a ^ 1 and the perfect non-polarity of the spacer, 
6 ^ 0. In addition, if we consider the dissociation energy 
of an oxygen from the surface Eg, which is the energy 
cost to generate the carrier at the surface, then the total 
energy is obtained as 

E = (L+toao+'71bo9)c^ — (i+2m/io)cx l+£'sC+const.. 

(4) 

Differentiating eq. (4) with respect to c, we obtain the 
optimum number of the carrier density as 

1„ L + 2mAo-Eg \ 
niAo + m_B02 / 

This optimum number c decreases from 1/2 (in the case 
without the spacer) because of the energy cost term of 
Eg. The other factors, the kinetic term, the hybridiza- 
tion effect at the interface, the lattice relaxation, and 
the atomic exchange, also decrease the carrier density at 
the interface. These effects will be discussed in the last 
section. 

2.2 Doping into both Interface and Spacer 

Next, to consider the effects of doping into the spacer 
in detail, we qualitatively estimate the optimum amount 
of the carrier in the substrate and spacer region again in 
the classical limit. We consider the case where the only 
one spacer is inserted above layer. The total energy 
of the interface mainly consists of 3 terms, the electric 
interaction term, the energy level of sites, and the dis- 
location energy. The electrostatic energies can be calcu- 
lated as a capacitor, where the energy is proportional to 
the thickness of the capacitor and the square of accumu- 
lated charge. We consider the hole doping into a p-type 
interface as well as the electron doping into the n-type 
interface (see Fig. 3). Then the total energy is given by 

E = {L + msos)^^ + Mbo2S^ + ("^B02 + Afsos)^^ 
-I- 2MBO2SC - 2(tobo2 + MB02)ch - 2MB02sh 
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Fig. 3. (Color online) Schematic of our model. Doped carriers are 
confined near the substrate (with the concentration — c) and near 
the spacer region (with the concentration —s). We also consider 
hole doping (with the concentration h) at the p-type interface 
between the spacer and the polar layer. 



-Lc-M B02 s+M BO2 h+Vs+V'h+Es (c+s- 



ft,)+const., 

(6) 

where MBO2 is the number of the polar layers above of 
the spacer, V and V' are the energy levels of the sites 
in the spacer region measured from the level of the sub- 
strate, s is the doped carrier density at the n-type in- 
terface of spacer, and h is the hole density at the p-type 
interface of spacer. To minimize the total energy eq. (6), 
the optimum carrier densities are obtained as 
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(9) 

where c, s, and h are > 0. We find that c and h in eqs. 
(7) and (9) do not depend on Eg. If Eg increases, the 
doped carrier s in the spacer region decreases first, and 
next, when s vanishes, the doped carrier c and hole h 
starts decreasing. A small change in Es has no effect on 
the carrier density at the interface. To consider the case 
/i = 0, a threshold V^j^ exists for the hole doing as 

(mB02 MBO2) 



th 



(10) 



Mbo2 + "^so2 

If V < V^^ is satisfied, the holes are not doped at the 
spacer. In this case, the optimum carrier densities are 
obtained as 

MBO2 + V 



1 L 
maxl-x — 



mB02 ~ Mbo2 



,0 



(11) 
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Fig. 4. (Color online) Schematic of our model. Filled circles rep- 
resent cations of non- polar substrate A^B^Os, open circles rep- 
resent cations of polar substrate B^Os and cross marked cir- 
cles represent cations of non-polar spacer A'^ B'^ O3. The doped 
electrons are assumed to transfer between the B sites in our 
model. 



and 



1 L~ Mbo, + V 
max — X 



2 L + mB02 - Mbo2 
1 , MBO2 - V -E^ 



2 2Mbo2 
In eq. (11), c does not depend on Es again. 



(12) 



3. Quantum Calculation by Hartree Approxima- 
tion 

3. 1 Formalism 

Next, we discuss quantum effects, namely, the effects 
of itinerancy of electrons, in our model. The carriers are 
expected to extend over several layers by the quantum 
effects. Moreover, if we consider the energy levels of the 
sites at the interface, the carriers are expected to be 
doped not only into the interface region, but also into 
the spacer region. 

Figure 4 illustrates a heterostructure which we con- 
sider in this paper. As shown in Fig. 4, B^O^ has 
a polarity, while the substrate and the spacer 

A^ B^Oy, do not. Each layer is a band insulator in the 
bulk. We neglect the orbital degrees of freedom of the 
transition metal d orbitals. We also neglect charge fluc- 
tuations of the valence bands of the oxygen 2p orbitals 
in this calculation. Our simplified Hamiltonian of this 
system is written as 

= Hhop + ^ ^pot + ^ surface? 

(13) 

i 

with i being the site index. The hopping term Tihop is 



defined as 



E 



(14) 



where cj^ (Qo-) is the creation (annihilation) operator of 
an electron on the ith-site with spin a. The summation 
over the nearest-neighbor sites is represented by {ij). The 
potential term Ti^j consists of three terms; 
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The one-body energy level is defined as 
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(15) 



(16) 



where rii = (^ia^ia is the number operator of the 
doped electrons. 

The long-range Coulomb interaction is defined as 
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e rijUi 
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(17) 



where the first and second terms are the Coulomb po- 
tential from the polar layer. Here, A^ and B^ represent 
the A sites and the B sites of the polar layers A^ B^O-^ 
respectively, and the third term comes from defects of 
oxygen. The last term is Coulomb repulsion between the 
doped electrons. For simplicity, the dielectric constants 
are assumed to be uniform in all the regions of the het- 
erostructure and the charge of the polarized layers are 
assumed to be ±e while an oxygen defect is assumed to 
generate the charge — 2e. The on-site Coulomb potential 
term is defined as 



on-site 



ViU. 



(18) 



The reconstruction energy of the surface is expressed by 

1 



surface — ^ ^ ^^-^si^'^^'^i 



(19) 



where Es{i) is the dissociation energy of an oxygen from 
the surface. In this study, we assume that the oxygens 
dissociate only from the single layer at the surface. From 
the electroneutrality condition, the carrier density and 
the oxygen-defect density satisfy the following condition: 



IE- 



Thus, eq. (19) can be written as 



rfacc 



= E. 



E 



(20) 



(21) 



where E^ has a unique value in each material. 

To consider the quantum effects quantitatively, we use 
the following parameters in our calculation. Lattice pa- 
rameter a is fixed at 3.9A, which corresponds to the ex- 
perimental lattice parameter of the SrTiOs substrate.* 
The hopping parameter t is set to 0.3 eV estimated from 
the bulk materials. ""^^ We take on-site Coulomb interac- 
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Fig. 5. (Color online) Total carrier density c as a function of 
dissociation energy Es- Total number of layers is taken as 50. 
Number of the polar layers L is indicated in the figure. By in- 
creasing the number of polar layers, the doped carrier densities 
increase. Increase in Eg leads to decrease in the carrier density 
from c = 0.5. 
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Fig. 6. (Color online) Carrier density n as a function of distance 
from surface d. Total number of layers is fixed at 50. Number of 
the polar layers L deposited between the surface and the interface 
is indicated in the figure as L/50. Dissociation energy Es is fixed 
at 154. The doped carriers are confined around the interface layer 
in the side of non-polar layers. 



tion in the substrate C/substratc = 6 eV (= 18t) estimated 
from high energy spectroscopies.^^ We employ the di- 
electric constant e = 15 and dimensionless parameter 



e /e ■ a ■ t 



0.8. The effect of the dielectric constant 



will be discussed in the last section. The energy level 
in polar layer Vpoiar is assumed to be 4.5 eV {— 15t) 
as is estimated from the bulk LaAlOa band structure. 
We neglect the on-site Coulomb interaction in the polar 
layer (?7poiar — eV) because the Al 3s wave function in 
LaAlOa/SrTiOs is widely extended compared to d and / 
orbitals. We leave the value of the dissociation energy Es 
as a parameter and examine the Eg dependence in the 
range ~ 5 eV. Different choices of Eg only affect the 
carrier density quantitatively and make no change in the 
essential mechanism of doping, if the doped carrier c is 
positive. The on-site Coulomb interaction in the spacer. 



Us, 



is also ignored for the sake of simplicity. To study 



the properties of this model, we employ the Hartree ap- 
proximation: Ui^n^i ^ (f^^T)"•ii - {nif){n,i). 
The results for the 100 x 100 x 50 site will be shown 
below. 




Fig. 7. (Color online) Total carrier density c vs. number of spacer 
m. Energy level Vspacer is fixed at I5t. Here, Es is assumed to be 
15t. Total number of layers is assumed to be 50 and the number 
of the polar layer is assumed to be 25. The spacers are inserted 
either between the BO2 and AO layers in the fifteenth polar 
layer or on the top of the fifteenth polar layer from the surface 
. Circle (red) plots show the result of the former case and open 
(green) squares show the latter. 



3.2 Quantum Calculation without Spacers 

We first calculate the charge distribution at the inter- 
face in the case without the non-polar spacer. The carrier 
densities of the interfaces are determined mainly by the 
balance between the potential energy and the dissocia- 
tion energy of the oxygen at the surface. As the thickness 
of the polar layers increases, the carrier densities tend to 
increase (see Fig. 5). This tendency is indeed found in 
real materials. Increase in the dissociation energy 
Es leads to decrease in the carrier density from c = 0.5. 

As shown in Fig. 6, the doped carriers tend to be con- 
fined near the interface. This tendency is also found in 
real materials. The polar instability is essentially de- 
termined by L, Es, e, and Vi, whereas the spread of the 
doped carrier extended away from the interface does nei- 
ther alter the instability nor the charge c. It only relaxes 



the local energy (we will discuss these effects in detail in 
the last section). 

3.3 Quantum Calculation with Spacer 

We now calculate the charge distribution at the inter- 
face in the case with spacers. The calculation was per- 
formed in a way similar to that in the previous subsec- 
tion. Figure 7 shows the doped carrier density as a func- 
tion of thickness of the spacer. The spacer is inserted 
either between the BO2 and AO layers in the fifteenth 
polar layer (contributing to mAo), or on the top of the 
fifteenth polar layer counted from the surface (contribut- 
ing to rn 302)- In the former case, the insertion of only 
ruAo = 10 spacers into 25 polar layers increase the carrier 
density up to around 50 percent compared to the non- 
inserted case. Furthermore, by the other type of spacers 
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Fig. 8. (Color online) Spacer-position dependence of carrier den- 



sity n. Here, Vsj 



is fixed at eV. Dissociation energy Eg is 



fixed at 15t. The number of total layers, polar layers and the 
thickness of spacer are taken to be 50, 25, and 5, respectively. 
The spacers are inserted on the top of BO2 layer of the polar 
layer. "Insert position" specifies the position of spacer counted 
from the surface. 



with TO_B02 1 the carrier density may not only be increased 
but also be decreased. In this way, the carrier density can 
be finely tuned similarly to the classical estimate in eq. 
(5). 

Next, we study the dependence of the doped carrier 
densities on the positions of the spacer inserted. Here, 
as an example, the spacer is inserted above the BO2 
layer, and the number of polar layers and the thickness 
of spacer are fixed at 25 and 5, respectively. The energy 
level T^paccr is fixed at eV. Figure 8 shows how the 
carrier density depends on the inserted position of the 
spacer. 

A characteristic feature is that, when the spacer is in- 
serted far from the surface, the carrier doping occurs in 
the spacer region as well. Whether the carriers are par- 
tially doped into the spacer region is determined by the 
position of the spacer. If the spacer is placed farther from 
the surface, the thickness of the polar layer between the 
spacer and the surface becomes larger and lowering of 
the energy by the carrier doping into the spacer region 
also becomes larger (see Fig. 8). As a result, carriers are 
doped partially into the spacer region. Of course, in the 
limit where the energy levels of the sites in the spacer 
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Fig. 9. Schematic of heterostructure for the case deviated from 
the ionic limit. Here, ±p and ±s represent charges of the layers. 
Carriers at the interface are doped not only from the surface of 
the polar layer (with the concentration — ci) but also from the 
bottom of the substrate (with the concentration — C2). 



region is sufficiently high, the carrier doping occurs only 
around the normal interface region. In this calculation, 
we neglect roles of the oxygen valence bands in the spacer 
because their energy levels are too low. 

In contrast to the electron doping, the lowering of en- 
ergy by the hole doping into the spacer region is enhanced 
when the spacer is not located near the interface but near 
the surface. This tendency is in contrast to the electron 
doping into the spacer. With this opposite tendency kept 
in mind, we should adjust the position of the spacer in 
order to control the carrier doping into the spacer. 

4. Why are Carrier Densities so low at the real 
Transition-Metal Oxide Interfaces? 

In this section, we examine the reasons why carrier 
densities are much lower than the expectation from the 
polar discontinuity in the available experimental condi- 
tions at the real transition-metal oxide interfaces. Later, 
we will also show that our doping method is consistent 
even in such low density cases seen in the experiments. 

Several works have found carrier densities of n-type 
interfaces around 10^'^ cm^^ much lower than the predic- 
tion considered in this paper, when the density of oxygen 
vacancies is expected to be low.^^ In fact, this is more 
than one order of magnitude less than 0.5 electron per 
unit cell derived in this paper. 

There are mainly three reasons for the low carrier den- 
sities at the interface. One of them is the insufficient 
thickness of the polar layers of the heterostructures. We 
have shown the thickness dependence of the doped carrier 
densities in the previous section. In the interface without 
the spacers, the doped carrier density c is obtained as 



max I i X (1 - ^),0 



(22) 



where e is the dielectric constant of the polar layers, Es 
is the energy cost of reconstruction and L is the number 
of the polar layers. If L is small compared to e and Es-, 
the carrier density c can be low or zero. This thickness 
dependence of carrier densities is negligible in the limit 
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of sufficiently thick L ^ oo. 

Even in such a thick limit, defects and substitutions of 
atoms may still have influences on the carrier densities 
at the interfaces. In the case where the interface has ad- 
ditional charges of electronic carriers c and atoms c', the 
total energy is obtained as 

E = L{c + cf - L{c + c') X 1 + const., (23) 

and the optimized carrier density is obtained as 

c = i - c'. (24) 

The instability of the potential divergence is removed by 
the reconstruction of both the electron and lattice sys- 
tems (see eq. (23)). If the most part of the reconstruction 
takes place in the lattice system, electronic carriers doped 
into the interface may become few. 

From the perspective of the polar discontinuity, the 
deviation from the ionic limit of the substrate also af- 
fects the expectation value of the doped carrier at the 
interface. In the ionic limit, the total energy is obtained 
in the classical electromagnetism as 

E = Lc^ - Lc + const., (25) 

where c is the doped carrier density per unit cell at the 
interface and L is the thickness of the polar layer along 
the [001] direction. Difl[erentiating eq. (25) with respect 
to the carrier density c, we obtain the optimum number 
of the carrier density as c = 1/2. On the other hand, if 
the system is away from the ionic limit, the total energy 
is obtained as 

E = Lcl+ L'cl - Lpci + L'sci + const., (26) 

where p is the charge of the polar bilayer (-|-p, — p), s is 
the charge of the "polarized" substrate bilayer (-|-s, — s), 
ci and C2 are the amounts of the charge at the surface 
of the polar layer and the substrate, respectively, and L' 
is the thickness of the substrate along the [001] direction 
(see Fig. (9)). The carrier density at the interface, c is 
obtained as the sum of ci and c^. Differentiating eq. (26) 
with respect to the carrier density c\ and C2, we obtain 
the optimum carrier density as 

c=\{3>-s). (27) 

We estimate charge of the LaO-A102 in the bulk LaAlOa 
and SrO-TiOa in the bulk SrTiOs by the local density 
approximation. The calculations were carried out with an 
in-house code, which is based on the Full-Potential Lin- 
ear Mufiin-Tin Orbitals (FP-LMTO) method.^" The lat- 
tice parameters are fixed at 3.905 A, which corresponds 
to the experimental lattice parameter of the bulk SrTiOa. 
In our calculations, LaO-A102 has the charge distribu- 
tion (-M, -1) and SrO-Ti02 has (-1-0.2, -0.2). Similar 
results of the charge distributions have been calculated 
by the generalized gradient approximation (GGA).^^ Us- 
ing these results, c is expected to be reduced by 20 per- 
cent from 1/2 even in the thick and clean limit. The 
strains from the substrates also affect the carrier densi- 
ties through the change in the ionicity of the polar layers. 
These three factors, the thickness, the lattice defects 
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and the deviation from the ionic limit, mainly affect the 
doped carrier densities at the interfaces. Other factors 
at the interfaces, atomic exchange, hybridization, carrier 
distribution and lattice relaxation, can also reduce total 
energy of the system and affect the doped carrier densi- 
ties.^^ However, in the limit where the thickness of the 
polar layers is sufficiently large (L oo), these local ef- 
fects do not alter the total accumulated charge around 
the interface. For example, the lattice relaxation at the 
interface changes the local polarity and the hybridiza- 
tion between the sites. However, the potential divergence 
along the [001] direction decreases only near the inter- 
face. The system in any case requires additional charges, 
the doped carriers or the defects, to avoid the instability 
of the potential energy. Quantitative explanation for the 
substantially small concentration of carriers observed ex- 
perimentally may be given by a combination of the three 
complex factors together with insufficient thickness of the 
LAO layers in the available experimental conditions. 

Our method of carrier control depends on the ratio 
of the thickness of the spacer to that of the polar layer. 
Even if the carrier densities decrease by the defects and 
the deviation from the ionic limit, the reduced carrier 
densities are still determined by the thickness of the polar 
layers and spacers, and our method can be used for the 
control of the carrier densities at the interfaces. 

5. Summary 

We have proposed a method to systematically control 
carrier densities at interfaces of transition-metal oxide 
heterostructures without introducing disorders. By in- 
serting non-polar layers sandwiched by polar layers, con- 
tinuous control of the carrier doping into the interface 
can be rc^alizcd. This method enables us to control the 
total carrier densities per unit cell c systematically up 
to high values of the order unity. We have demonstrated 
this charge controllability and the resultant distribution 
of the confined carriers at the interface by using classical 
and itinerant models in the Hartree mean-field theory. 
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